Bridge relations in Navier-Stokes turbulence 



00 



> ■ 

i> : 

(N . 
o ■ 
\o • 

o : 

00 . 

6 : 



A. Celani ^-^ and D. Biskamp ^ 
^ Max Planck Institut fiir Plasmaphysik, D-8574-8 Garching, Germany 
INFM Unita di Torino Universitd, Italy 

The correlation between inertial range velocity fluctuations and energy dissipation in fully devel- 
oped turbulence is studied using high resolution direct numerical simulation. Runs with microscale 
Reynolds number up to TZ\ ~ 300 have been performed in order to achieve the maximum extent of 
the scaling range compatible with a proper resolution of the dissipative scales. The bridge relations, 
a set of scaling relations which express quantitatively the statistical dependence of dissipation-scale 
velocity fluctuations on inertial-range ones, are verified. Corrections to the bridge relations due to 
the finite size of the scaling range are investigated. 
PACS: 47.27.Gs, 47.27.Eq 



Fully developed Navier-Stokes turbulence is character- 
ized by an energy flux from large scales, L, driven by 
some stirring force, to small scales, 77, where viscous dis- 
sipation is taking place. Within the intermediate range 
of scales L 3> r ^ 77, called the inertial range, in which 
neither forcing nor dissipation are active, longitudinal ve- 
locity differences Sru{x) = {u{x + r) — u{x))-r /r display 
a scaling behavior: in homogeneous and isotropic turbu- 
lence, structure functions 



Sp{r) - {{6ru{x))P) 
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behave as power-laws Sp{r) ~ r^p. According to the 
Kolmogorov 1941 theory = p/3, on the basis of di- 
mensional arguments and self-similarity Actually 
there is strong evidence that exponents Cp are not those 
predicted by Kolmogorov, but are a nonlinear function 
of p (see e.g. |^). This deviation from the Kolmogorov 
1941 prediction is the most important signature of inter- 
mittency in the energy transfer process. 

A basic consequence arising from the existence of an 
energy transfer directed to small scales is the statistical 
dependence of small-scale velocity fluctuations on large- 
scale ones. This dependence can be quantitatively in- 
vestigated by looking at the behavior of inertial range 
two-scale velocity correlations such as 



(2) 



with L R r rj. Under the usual scaling hypothe- 
sis and a mild universality assumption on the small-scale 
structure of turbulence, it can be shown that, to leading 
order, |,| 



Fp^,{r,R)^Sp+,iR) 
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with an exception for p = 1, where the prediction (|^) 
does not apply due to the fact that Si{r) vanishes be- 
cause of homogeneity (see Refs. These relations, 
and similar ones that have been derived in the same con- 
text, have been successfully tested in experiments 
Since the scaling relations (||) hold exactly only when 



£ 3> i? 3> r ^ ?7, their assessment in numerical simula- 
tions is currently out of reach due to the limited extent of 
the inertial range achieved with the actual computational 
capabilities. It is worth noting that the result (||) agrees 
with the prediction of a multiplicative random process 
for the energy transfer mechanism 0| . 

Since the energy flux towards small scales will flnally 
be dissipated by viscous shear, it is reasonable to ex- 
pect that, since this transfer induces correlations be- 
tween inertial-range fluctuations at different scales, it will 
also induce correlations between inertial-range fluctua- 
tions and dissipative-scale ones. This is the same kind 
of physical picture that lies at the base of Kolmogorov's 
Refined Similarity Hypothesis 

A recent result is that quantitative predictions about 
the correlation between energy dissipation s{x) = 
z^l VM(a;)p and inertial range velocity fluctuations Sru{x) 
can be derived from Navier-Stokes equations under the 
same hypotheses used in the derivation of (||) These 
quantitative predictions involve correlations of the kind 

D^'Hr) ^ {e{x)i6Mx)r) , (4) 
D(2)(r) = {£{x)e{x + r)(Sru{x))P) , (5) 

which have been shown to follow the bridge relations M 
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The hierarchy of bridge relations extends to higher or- 
der correlations involving three or more factors with 
energy dissipation. Let us stress that in the presence 
of intermittency the scaling relations (Hj^) are differ- 
ent from the scaling derived if dissipative scale fluctu- 
ations were independent of inertial scale ones. Indeed 
the independence of e{x) from SrU would have led to 
Dl^\r) ~ (e)S'p(r) ~ r^p; this is in contrast to (^), from 
which follows that D'j^\r) 



■Cp + 3 



smce m presence 
of intermittency we have that Cp+3 ~ 1 < Cp- 

Note that the nature of scaling relations (0) - which are 
derived only under general assumptions on scaling and 
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universality of small-scale turbulence, and are shared by 
all systems displaying a direct energy cascade - is very 
different from the nature of the bridge relations 
which are derived from the Navier-Stokes equation and 
thus have a dynamical content. It is noteworthy that the 
bridge relations coincide with the predictions given by 
a multiplicative random process with a proper account- 
ing of the viscous cutoff scales, and making explicitly 
use of the Kolmogorov 4/5 law in the form Cs = 1 0- 
The actual need for a relation like Kolmogorov's law, 
which is itself a consequence of the Navier-Stokes equa- 
tion, to derive (^|j^) in the framework of multiplicative 
models shows again the inherently dynamical nature of 
the bridge relations. 

Among the bridge relations there are two renowned 
special cases. The first one is ^ for p = 0, which gives 
(e) - S3{r)/r. This is the Kolmogorov 4/5 law 
which is the only bridge relation that can be established 
within coefficients 
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The second well-known bridge relation is (|^) for p = 0, 
{e{x)e{x + r)) ^ Se{r)/r'^. It implies the scaling for the 
two-point dissipation correlation {e{x)e{x + r)) ~ r~'^ 
with /Li = 2 — (^6- This relation involving the intermit- 
tency exponent fi is the phenomenological "bridge rela- 
tion" |^,||, hence the name chosen for (^,0). 

In this Letter we present the results of direct numeri- 
cal simulations of fully developed homogeneous isotropic 
turbulence which show clear evidence for the validity of 
the bridge relations (^,0) in the inertial range of scales. 

The simulations solve the Navier-Stokes equation on 
a grid with 512^ collocation points in a periodic cubic 
box of linear size 27r, using a standard pseudo-spectral 
method. Time stepping is performed with a trapezoidal 
leap-frog scheme and the dissipative term is integrated 
exactly. To obtain a statistically stationary state the ve- 
locity field is forced for wavenumbers k < 2.5. The forc- 
ing scheme keeps the total energy constant in each of the 
first two shells (0.5 < k < 1.5 and 1.5 < k < 2.5) with a 
ratio consistent with k^^^^ jl^. Two runs with different 
viscosity have been performed, with microscale Reynolds 
numbers TZ\ ~ 300 and TZ\ ~ 220, defined in the usual 
way as TZ\ = vqX/v, where A = (IbiyvQ/eY^^ is the Tay- 
lor microscale and vq is the root mean square value of a 
single velocity component. With the Kolmogorov length 
r] — (i/^/e)^/^, the two runs correspond, respectively, 
to kmaxV ^ 1 and kmaxV ~ 2, where kmax = 256 is 
the largest wavenumber resolved. The run at the larger 
Reynolds number, which has been chosen to maximize 
the scaling range, resolves the dissipative scales only 
marginally well, a fact that might shed some doubt on 
the relevance of the results for the bridge relations (^U^) , 
since they involve the energy dissipation e{x). We have 
therefore also performed the lower Reynolds number run, 



which has adequate resolution of the dissipative range, 
and thus serves as a benchmark to test the quality of 
the data obtained at TZ\ ~ 300. Ensemble averages are 
replaced by spatial averages and time averages were per- 
formed over ten frames for each run in a time-span of 
approximately 7 large eddy turnover times. 

In Figure |l| is shown the second-order longitudinal 
structure function 52 (r) rescaled with respect to the 
square of the Kolmogorov velocity = e^^^rj^^^, for the 
two runs at different Reynolds number. 
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FIG. 1. Structure functions S2(r)/v^ for longitudinal ve- 
locity increments at different Reynolds numbers, boxes (□) 
for TZx ~ 300, triangles (A) for Tlx ~ 220. The solid line is 
the Batchelor fit ^ with b = 10.17. In the inset: the local 
slope C2(r) for the data at Tlx ^ 300 (□) and Tlx ^ 220 (A), 
compared with the value extracted by ESS ("2 ~ 0.69 (solid 
line) and the Kolmogorov value ("2 = 2/3 (dashed line) 

According to the Batchelor parametrization, 5*2 (r) can 
be expressed as |p|,|ll| 



S2ir) 



[l + (l + l/36)3/2(r/r/)2]i-C2/2 



(9) 



which indeed fits the data very well, when the value 
C2 = 0.69 extracted by ESS technique is assumed [T^ . 
The almost perfect agreement of the data for TZ\ ~ 220 
and TZ\ ~ 300 at small and intermediate scales demon- 
strates the relevance of results of the run at larger 
Reynolds number. In the inset of Figure |^ is also 
shown the local slope of the second order structure func- 
tion C,2{'r) ~ d\nS2{r)/d\nr. Notwithstanding the lim- 
ited scaling range, for TZ\ ~ 300 an interval in which 
C,2{'r) is approximately constant can be detected, with 
C2(?') — 0.69, definitely larger than the Kolmogorov value 
2/3 and in agreement with the value obtained using ESS. 
In the case with TZ\ ~ 220 no definite slope can be in- 
ferred from C2(f), and one has to rely completely upon 
ESS to determine the scaling exponent. 
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In Figure || we show the longitudinal structure func- 
tions for the lowest even orders, and a check of Kol- 
mogorov 4/5 law (|), both for Tlx ^ 300. 
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FIG. 2. Structure functions Sp{r) for longitudinal velocity 
increments with p = 2,4, 6. Solid lines are power laws with 
exponents C2 = 0.69, C4 = 1-27,(6 = 1.77 extracted by ESS 
technique. Data are shifted on the vertical axis for viewing 
purposes (factor lO'^ for p = 4 and 10'' for p — 6). In the inset 
is shown the Kolmogorov 4/5 law (^). Note the linear scale 
on the vertical axis. 

These results are in agreement with previous numerical 
simulations p^ ]. 

In Figure ^ we examine the scaling behavior of the 
correlations DjP{r) — {e{x){dru{x))P). 
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FIG. 3. Correlations Dp'' (r) for longitudinal velocity incre- 
ments with p = 2, 4. Solid lines are power laws with exponents 
(5 — 1 = 0.53, (7 — 1 = 0.99, which are the scaling exponents 
predicted by the bridge relations (^). The exponents (p have 
been obtained by ESS technique. 



According to the bridge relations (H), the scaling 
D^\r) ^ should hold within the inertial range. 



The data are in good agreement with the theoretical pre- 
diction within the accuracy allowed by the limited scaling 
range. 

In order to have a more accurate assessment of the 
bridge relations (||), it is useful to display the behavior 
of the ratio rDp^\r)/Sp+3{r), which should be indepen- 
dent of r within the inertial range. As shown in Figure |^ 
there is a clear plateau (note the linear scale on the ver- 
tical axis) in an interval of values of r coincident with 
the scaling range of the even-order structure functions, 
confirming the results of Figure ^. 
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FIG. 4. The ratio rDl^\r)/Sp+:i{r), approximately inde- 
pendent of r for intermediate values of r, verifying the bridge 
relations (^). Odd order Sp{r) are computed taking the ab- 
solute value of velocity differences. 

Figure ^ also shows that the bridge relations break down 
at small and at large r. Indeed close to and below 
the dissipative scale 77 we have the smooth behavior 



Dl^>{r) ~ (eP/2+i)rP/i.P/2 and ^p(r) 



.p/2 



which leads to the scaling relation 

D^^\r) ^ i^Sp+2{r)/r^ for r/rj . (10) 

Furthermore, when the distance r becomes comparable 
with the integral scale L we expect that the statistics of 
velocity differences becomes independent of the dissipa- 
tion. Thus the correlation d'^^ (r) must asymptotically 
become 



Di'Hr)^{e)Spir) 



for 



r > L 



(11) 



The asymptotic behavior for Dp^"* (r) is shown in the up- 
per panels of Figure ^. 

In order to check the validity of the bridge relations M) 
we plotted the ratio r^Dp {r) / Sp+f,{r) versus r, which - 
in agreement with (|^) - is constant for values of r within 
the inertial range of scales, as shown in Figure I The 
case p = corresponds to the phenomenological bridge 
relation (e(a;)e(a; -I- r)) ^ ~ r^"?^, with /i = 0.22. 
The agreement is slightly less good than that for D^'' (r) , 
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an effect that can be ascribed to the higher order of the 
structure functions involved. 
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FIG. 5. The ratio r'^D^\r)/Sp+Q{r), approximately inde- 
pendent of r for intermediate values of r, verifying the bridge 
relations (^). 

For Dp (r) there are also departures from the bridge re- 
lations approaching either the dissipative scale rj or the 
integral scale L, with leading scaling laws 

D^^\r)r.v''Sp+i{r)/r^ for r/77 ^ (12) 

Df\r)^{e)^Sp{r) for r>L, (13) 

which follow respectively from the smoothness of the ve- 
locity field at dissipative scales and from the indepen- 
dence of the dissipative range on the forcing scales, and 
are shown in the lower panels of Figure ^. 
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FIG. 6. Upper Left: The ratio r'^ D^p\r) / {vSp+2{r)) be- 
comes independent of r for r/77 0, according to (p^. 
Upper Right: The asymptotic decorrelation of energy dis- 
sipation and velocity differences reflects in the behavior of 
D^\r) / {£Sp{r)) which tends to unity for r/L ~ 1. Lower 
Left: r"^ D'^^t) / {u^ Sp+4,{r)) becomes independent of r for 



77, according to 



Lower Right: Asymptotic behavior 



In summary, the physical picture emerging is that the 
straining process transforms inertial scale fluctuations 
into dissipation scale fluctuations, thus inducing a statis- 
tical dependence of the latter on the former with an in- 
herently dynamical nature. The "bridge relations" are a 
quantitative expression of this dependence, which deter- 
mine the scaling exponents of correlations between dissi- 
pation and inertial range velocity fluctuations in terms of 
the usual scaling exponents of structure functions. These 
relations have been successfully tested in the framework 
of direct numerical simulations of Navier-Stokes turbu- 
lence. The analysis has been performed with microscale 
Reynolds number up to TZx ~ 300, a value that allows to 
achieve the maximum extent of the scaling range without 
giving up the proper resolution of the dissipative scales. 
Experimental confirmation of the results here presented 
would be heartily welcome. 
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